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Abstract. We calculate improved lower bounds for the connective constants for self- 
avoiding walks on the square, hexagonal, triangular, (4.8 2 ), and (3.12 2 ) lattices. The 
bound is found by Kesten's method of irreducible bridges. This involves using transfer- 
matrix techniques to exactly enumerate the number of bridges of a given span to very 
many steps. Upper bounds are obtained from recent exact enumeration data for the 
number of self-avoiding walks and compared to current best available upper bounds 
from other methods. 
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1. Introduction 

The self-avoiding walk (SAW) on regular lattices is one of the most important and 
classic combinatorial problems in statistical mechanics ^Hl- An n-step self-avoiding 
walk w on a regular lattice is a sequence of distinct vertices uj ,ui, . . . ,u n such that 
each vertex is a nearest neighbour of it predecessor. SAWs are considered distinct up to 
translations of the starting point Uq. The fundamental problem is the calculation (up 
to translation) of the number of SAWs, c n , with n steps. It is generally believed that c n 
grows exponentially with power law corrections 

Cn ~ vw- 1 , (1) 

where \i is called the connective constant, 7 is a critical exponent and A a critical 
amplitude. Hammersley and Morton JU] were the first to prove the existence of the 
limit 

/1 = lim <#» (2) 

n— »oo 

The exact value of fi is known only on the hexagonal lattice, where Nienhuis [T§1 |2*U] 
showed, using non-rigorous methods, that //hex = v2 + y/2, and on the (3.12 2 ) lattice, 
where Jensen and Guttmann ^S] found an exact and rigorous connection between 
the connective constant /i(3.i2 2 ) and the connective constant for the hexagonal lattice 
A^hex = A^ 3 i 2 2)/(A i (3.i2 2 ) + 1)- Cm the square lattice it has been observed [H] that /x sq 
is indistinguishable from the reciprocal of the unique positive root x c of the simple 
polynomial 581x 4 + 7x 2 — 13 = 0, and while this 'conjecture' has stood the test of time 
it remains a purely numerical observation. 

Since finding the exact value of \i (let alone proving such results rigorously) is 
extremely difficult much effort has been devoted to more general methods for proving 
rigorous bounds on the connective constant. Brief overviews of some of the methods 
used can be found in [TBI E! • A systematic procedure for improving the lower bounds 
can be devised from a method due to Kesten [TH]. It was used by Guttmann [7j to 
improve the lower bounds for the connective constant on the square and simple cubic 
lattices and more recently by Aim and Parviainen [3] to obtain improved lower bounds 
on the connective constant for the hexagonal lattice. In this paper we further refine 
these bounds and extend the work to the triangular, kagome and (4.8 2 ) lattices. 

Finally, we use recent exact enumeration data for c n to obtain upper bounds for 
the connective constant on the square, hexagonal and triangular lattices. These bounds 
are then compared to better upper bounds obtained from other methods. 

2. Lower bounds 

Lower bounds for the connective constant can be found using a method due to Kesten 
[TT)] . The method utilises the fundamental result that certain restricted classes of 
self-avoiding walks have the same connective constant as the unrestricted problem. 
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Particularly useful for our purposes is the class of walks known as bridges. Let Xj denote 
the x-coordinate of Uj, then a bridge is a self- avoiding walk such that xq < Xj < x n for 
all j > 0. We use b n to denote the number of n-step bridges, and note Kesten showed 
that b\[ n converges to fi as n — > oo. Clearly concatenating two bridges of length n 
and m gives a bridge of length n + m (we place the origin of the second walk on top 
of the end-point of the first walk). This means that any bridge can be decomposed 
into irreducible bridges, i.e., bridges which cannot be decomposed further, and we use 
a n to denote the number of n-step irreducible bridges. It is now easy to see that the 
generating function B(x) for bridges is simply related to the generating function for 
irreducible bridges A(x) 

B{x) ' 



1 - A(x) 

It then follows that 1/fi is the solution to A(x) = 1. This relation also allows us to 
obtain lower bounds for \i. This relies on the observation that, if < a n < a n , for 
n > 2, then with x c being the solution to 

oo 

J2*nX n = l (3) 

n=l 

l/x c is a lower bound on \i. In particular we can set a n = for n > N and thus truncate 
the series. 

It is not easy to calculate the number of irreducible bridges directly. Thankfully 
they can easily be obtained from the number of bridges. Following Aim and Parviainen 
|lj we consider the number of bridges b n j and irreducible bridges a n j of length n and 
span /, that is bridges with xq = and x n = I > 0, with associated generating functions 
B[(x) and Ai(x). Obviously YlHi a n,i — a n, so if we truncate at some maximum span L 
and maximum walk length iV then the reciprocal of the solution to 

N L 

j2J2 a ^ n = 1 ( 4 ) 

n=l 1=1 

is a lower bound on 

Since a bridge is either irreducible or the concatenation of a bridge with an 
irreducible bridge we get 

i-i 

Bi(x) = Mx) + J2 A i-k(x)B k (x) 

and thus 



k=l 



i-i 



Ai(x) = Bi(x) -J2 A i-k( x )B k (x), 

k=l 

which allows to obtain all generating functions A\{x) recursively from Bi(x) for 1 < I < 
L. 
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In this paper we also examine a second way of obtaining lower bounds. We again 
use irreducible bridges, but rather than using smallish L and very large N we calculate 
the exact series for irreducible bridges to order N (much lower than before) and use this 
truncated series to obtain a lower bound from the reciprocal of the solution to 



that is in Eq. (JSJ) we set a n = a n for n < N and a n = for n > N. 
2.1. Enumeration of self- avoiding bridges 

The number of self-avoiding bridges b n j can easily be counted using the Transfer-Matrix 
(TM) methods we have developed for the unrestricted problems O EH E], which 
are devised to count the number of walks in a finite I x w rectangular sub-section 
of the underlying lattice. Here we shall only briefly outline the changes required to 
enumerate bridges. The most efficient implementation of the TM algorithm generally 
involves bisecting the rectangle with a boundary line and moving the boundary in such 
a way as to build up the lattice cell by cell. The sum over all contributing graphs 
is calculated as the boundary is moved through the lattice. For each configuration of 
occupied or empty edges along the intersection we maintain a generating function for 
partial walks cutting the intersection in that particular pattern. If we draw a SAW and 
then cut it by a line we observe that the partial SAW to the left of this line consists 
of a number of loops connecting two edges in the intersection, and at most two pieces 
connected to only one edge (these are the pieces from the end-points of the SAW). The 
computational complexity of the algorithm is essentially determined by the number of 
such configurations. So we must make the intersection as short as possible. Since we 
are looking to fix I and make N large it follows that w will be large as well (in fact 
proportional to N — I). So the boundary line must intersect the rectangle along the 
'bridging' axis, e.g., along up to I + 2 edges. It is quite easy to demonstrate 0112] that 
the number of configurations grows like 3' in the square lattice case. So the required 
CPU time will grow roughly as (w + l) 2 3 l = N 2 3 l , since there are (w + 1) — N updates 
and terms in the generating functions. Memory requirement will grow as N3 l . We note 
in passing that while the TM method can be used to study higher-dimensional lattices 
it quickly becomes inefficient because the boundary would be (d — 1) dimensional and 
the number of edges in the intersection would grow ever more rapidly. 

In order to implement the first method for finding lower bounds, see Eq. (j3J), we 
count the number of bridges spanning rectangles of size I xw, that is bridges starting at 
the bottom border and terminating at the top border. In addition the walks must also 
touch the left border of the rectangle (this takes care of the translational invariance) as 
illustrated in Figure H in all case bridges must terminate at a top most vertex in the 
top most row. Note in particular the implication of this restriction on the hexagonal, 
kagome, and (4.8 2 ) lattices. For the hexagonal case this means that all bridges are of 
even length. 



iV 




(5) 



n=l 
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Figure 1. Examples of bridges on the square, hexagonal, triangular, kagome and 
(4.8 2 ) lattices. In addition we also show a section of the (3.12 2 ) lattice. 



For the square lattice we calculated the number of bridges to L = 15 and iV = 250, 
for the hexagonal lattice to L = 15 and N = 500, for the triangular lattice to L = 12 
and iV = 150, for the kagome lattice to L = 10 and N = 300, and for the (4.8 2 ) lattice to 
L = 12 and N = 500. Because of the exact connection between the connective constants 
A*(3.12 2 ) an d /Uhex, /^hex = A*( 3 i 2 2)/(A t (3.i2 2 ) + 1), an y bounds for the hexagonal lattice yields 
corresponding bounds for the (3.12 2 ) lattice. So we don't actually count bridges on the 
(3.12 2 ) lattice and have thus not shown an example of one in Figure ^ 

The integer coefficients occurring in the series expansions become very large. The 
calculations were therefore performed using modular arithmetic [T7]. This involves 
performing the calculation modulo various integers p, and then reconstructing the full 
integer coefficients at the end. The Pi are called moduli and must be chosen so they are 
mutually prime. The Chinese remainder theorem ensures that any integer has a unique 
representation in terms of residues. If the largest value occurring in the final expansion 
is m, then we have to use a number of moduli k such that P1P2 ■ ■ - pk > rri- We used 
moduli which are prime numbers of the form pk = 2 30 — r*.. 
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Naturally the calculation for each I and moduli are completely independent. It is 
evident from the exponential growth in the computational complexity most of the CPU 
time is spent on the largest value of L, where up to 16 moduli were required to represent 
the coefficients. Typically each run (at the maximal span L) required up to 24 CPU 
hours on a 1GHz Alpha processor and could use up to 2.5Gb of memory. In all we used 
about 3000 CPU hours on the calculations. Our method is much more efficient than 
that used by Aim and Parviainen 0], who report using more than 20000 CPU hours 
calculating the number of bridges on the hexagonal lattice with L = 10 and N = 58. A 
similar calculation using our method takes no more than a couple of minutes! 

The second method for finding lower bounds, see Eq. ( |5|) . uses the exact data for 
the number of irreducible bridges up to length N. Again the first step is the calculation 
of the relevant data for bridges. We illustrate the method in the square lattice case. An 
irreducible bridge of width L has length at least 3L. This is because each row (apart 
from the bottom most) must have more than one occupied edge (otherwise we could cut 
the walk into two bridges) and the walk must thus go up, come down, and go up again. 
We also have the first step and at least two horizontal steps for a grand total of at least 
3L steps. So if we require the number of irreducible bridges to order N we must count 
the number of bridges with span up to L = N/3. That is we have to count the number 
of bridges on rectangles of size w x I, where 1 < I < L — N/3 and 1 < w < N — I. Note 
that this calculation gives the number of bridges correctly only to order L. However, by 
first extracting the series for A(x) we can also get B(x) = 1/(1 — A(x)) correct to order 
3L. The efficient calculation of the bridge generating function is in many aspects more 
complicated and time consuming than for the first method. Details of the properties 
of the bridge generating function will appear in a separate paper. Suffice to say that 
we have obtained generating functions to order 72 on the square lattice and 122 on the 
hexagonal lattice. 

The series for the problems studied in this paper can be obtained by request from 
the author or at http: / /www.ms. unimelb.edu.au/~iwan/ by following the relevant links. 



2.2. Results 

Lower bounds are obtained by forming the polynomials of Eq. (j3J). It is thus possible 
to obtain ever improved lower bounds by increasing N and L. In Table ^ we use the 
hexagonal lattice data to illustrate the method. Note that for this problem the minimal 
number of steps for an irreducible bridge of span L is 6L — 2. From this data we 
observe first of all, that for N = 100 little is gained by going beyond span L — 12. 
This is somewhat surprising since A 13 (x) contributes already at order 76, but obviously 
the influence of these terms is almost negligible. Likewise with fixed L and increasing 
N it is a case of rapidly diminishing returns. If we are interested in optimising the 
procedure, that is, getting a decent bound, but with as little wasted effort possible, it 
appears that for given L we should choose N larger than twice the order of the first 
non-zero contribution to A L (x) (otherwise the calculation of A L (x) is largely wasted) 
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Table 1. Lower bounds for the conn ective constant for the hexagonal lattice. The 
supposed exact value is [i = y/2 + \/2 = 1.847759065 . . .. 



L N = 100 JV = 200 AT = 300 N = 400 = 500 



2 


1 


.787760 


1 


.787760 


1 


.787760 


1 


.787760 


1 


.787760 


3 


1 


.808678 


1 


.808678 


1 


.808678 


1 


.808678 


1 


.808678 


4 


1 


.819369 


1 


.819370 


1 


.819370 


1 


.819370 


1 


.819370 


5 


1 


.825750 


1 


.825767 


1 


.825767 


1 


.825767 


1 


.825767 


G 


1 


.829869 


1 


.829984 


1 


.829984 


1 


.829984 


1 


.829984 


7 


1 


.832546 


1 


.832950 


1 


.832951 


1 


.832951 


1 


.832951 


8 


1 


.834182 


1 


.835137 


1 


.835140 


1 


.835140 


1 


.835140 


9 


1 


.835061 


1 


.836802 


1 


.836814 


1 


.836814 


1 


.836814 


10 


1 


.835448 


1 


.838094 


1 


.838132 


1 


.838132 


1 


.838132 


11 


1 


.835574 


1 


.839101 


1 


.839191 


1 


.839193 


1 


.839193 


12 


1 


.835602 


1 


.839878 


1 


.840058 


1 


.840063 


1 


.840063 


13 


1 


.835606 


1 


.840464 


1 


.840775 


1 


.840788 


1 


.840789 


14 


1 


.835606 


1 


.840890 


1 


.841372 


1 


.841400 


1 


.841402 


15 


1 


.835606 


1 


.841184 


1 


.841868 


1 


.841921 


1 


.841925 



but not much larger than four times this order. Similar considerations apply to the 
other problems as well though the optimal cut-off varies from problem to problem. 

Here we briefly summarise our results for the lower bounds. For the hexagonal 
lattice we find the lower bound, 1.841925 < /ihex, which is less than 0.32% lower 
than the exact value /ihex = v2 + a/2 = 1.847759065 .... The previous best lower 
bound was 1.833009 < fi^ex 111- For the square lattice we obtain the lower bound, 
2.625622 < /i sq , which is within 0.48% of the best estimate for the connective constant 
fi sq = 2.63815853031(3) [TT]. This should be compared to the previous bound 2.62006 < 
/4q [6J- For the triangular lattice the lower bound is, 4.118935 < /x tr i, within 0.77% of 
the best estimate /i tr i = 4.150797226(26) [TH], whereas the previous best bound was 
4.03333 < 

/Vi 0- The Kagome lattice lower bound is, 2.548497 < /Xkag> within 0.48% 
of the estimate /ik ag = 2.560576765(10) (based on our unpublished enumerations of self- 
avoiding polygons), while the previous best bound was 2.50967 < ;U kag |3|. For the (4.8 2 ) 
lattice we found the lower bound, 1.804596 < /i(4.s 2 ), which is just 0.24% lower that the 
estimate — 1-80883001(6) [T3], which improves on the bound 1.78564 < /i(4.s 2 ) 0- 

Finally, for the (3.12 2 ) lattice we get the lower bound, 1.708553 < /i(3.i2 2 )) which is just 
0.15% from the exact value /i(3.i2 2 ) = 1.711041296. . . [THJ , and again improves on the 
previous bound 1.705263 < /i(3.i2 2 ) 111- 

As stated earlier we also used a second approach to obtain lower bounds for the 
connective constant. This entails the calculation of an exact series expansion for the 
generating function for irreducible bridges up to some maximal order N (this was also 
the method employed by Guttmann [HE])- Lower bounds are then obtained from the 
truncated series in Eq. (jSJ). Obviously we could truncate at any order n < N and obtain 
a sequence of lower bounds n{n). In Table El we have listed the lower bounds obtained 
from this method for the hexagonal and square lattice cases. Clearly this method 
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Table 2. Lower bounds for the connective constant for the hexagonal and square 
lattices. 



Hexagonal Square 



1 V 


1— </"n lvirl 

L>U LI 11 LI 


M 

1 V 


1— </~n mn 
L>U Li 11 Li 


32 


1.812833 


30 


2.583704 


38 


1.817977 


33 


2.588448 


44 


1.821786 


36 


2.592419 


50 


1.824722 


39 


2.595794 


ob 


l.oz /Uoo 


4z 


Z.o9oo9o 


62 


1.828955 


45 


2.601224 


G8 


1.830532 


48 


2.603442 


74 


1.831863 


51 


2.605405 


80 


1.833002 


54 


2.607155 


8G 


1.833987 


57 


2.608726 


92 


1.834847 


60 


2.610143 


98 


1.835606 


63 


2.611428 


104 


1.836279 


66 


2.612599 


110 


1.836882 


69 


2.613671 


116 


1.837424 


72 


2.614656 



is inferior to the previous one (the bounds are not as good) particularly considering 
that the computational effort is significantly greater. However, this approach allows 
us to study the convergence of the lower bounds fi(n) to the connective constant as a 
function of the truncation order n. We find that fi — fi(n) ~ a/n, this behaviour can 
be seen directly in Fig. |21 where we have plotted fi — fi(n) vs. 1/n. We also formed 
the generating function D(x) = ^2 n d n x n , where d n = /i — fi(n), analysed this using 
differential approximants and found a logarithmic singularity at expected if 

d n ~ a/n. 

3. Upper bounds 

The best current method for obtaining upper bounds is due to Aim [2] and it essentially 
requires one to enumerate the number of walks according to length n and a specified 
head and tail each of length m. More precisely Aim showed that 

At < (A(G(m,n))) 1/( "" m) , (6) 

where A is the largest eigenvalue of the matrix G(m,n). The entries of this matrix 
are equal to the number of n-step self-avoiding walks starting with a walk u)i and ending 
with a translation of a walk ojj. Each walk oj^ i — 1, . . . , K m , is one of the K m possible 
m-step self-avoiding walks (up to all possible symmetries). While this method can yield 
quite sharp upper bounds (within 1.1% for the hexagonal lattice jl]) it is unfortunately 
not suited for a transfer-matrix enumeration. 
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Figure 2. The difference between the connective constant /i and the lower bound 
/Lt(n) vs. 1/n. Data is plotted for the hexagonal and square lattices. 



However, we have recently obtained greatly extended series for the number of SAWs 
on the square, hexagonal and triangular lattices [T2j [TBI El- This allows us to use a 
special case of Aim's work j2] which states that if K m = 1 then 



fl< fl % 



n 



ipn/ Cm) 



l/(n—m) 



(7) 



On all lattices K\ = 1 so fi < fi\{n) = (c n /ci) 1 ^ n_1 - ) as proven earlier by Ahlberg 
and Janson pQ while K2 = 1 for the hexagonal and L lattices yielding sharper bounds 
fi < fi2(n) = (c n /c 2 ) 1 ^ n_2 - ) as proven by Guttmann |S]. Ahlberg and Janson also proved 
that an upper bound fi a {n) can be obtained from the positive root of the equation 



zx 



_n-l 



[c n -{z- 2)c n - 1 ]x + (z - 2)[(z - l)c n _i - c n ), 



(8) 



where z = c\ is the coordination number of the lattice. An upper bound is then found 
as min(/! m (n),/i a (n)). 

For the square lattice we have c 71 = 4190893020903935054619120005916 and c 70 = 
1580784678250571882017480243636, which gives us the upper bounds fi x = 2.684484 
and fi a = 2.681360. These bounds are sharper than the one obtained by Aim |2j, 
fi sq < 2.695759, using n = 24 and m = 8. An improved upper bound, fi sq < 2.679193, 
has been obtained by Ponitz and Tillmann [2JJ, by counting walks with finite memory. 

For the triangular lattice we have c 40 = 22610911672575426510653226 and c 39 = 
5401678666643658402327390, which gives us the upper bounds fi x = 4.267349 and fi a = 
4.263713. These bounds are sharper than the one obtained by Aim [2], fi t n < 4.277799, 
using n = 16 and m = 6 (Aim has since improved this to /x t ri < 4.25152 [3J. 
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Table 3. Summary of results for the connective constant [i with the current best 
lower and upper bounds 



Lattice 


Lower Bound 




Upper Bound 


Square 


2.625622 


2.63815853031(3) [TT| 


2.679193 [H] 


Hexagonal 


1.841925 


1.847759065... |H] 


1.868832 gj 


Triangular 


4.118935 


4.150797226(26) [B| 


4.25152 


Kagome 


2.548497 


2.560576765(10) 


2.590301 


(4.8 2 ) 


1.804596 


1.80883001(6) [IS] 


1.82926 


(3.12 2 ) 


1.708553 


1.711041296... [El 


1.719254 gj 



For the hexagonal lattice we have c 100 = 2585241775338665938539885252 and 
c 99 = 1394474897269109512317080364, which gives us the upper bounds ^ = 1.871004, 
H2 = 1.869731 and /i a = 1.869836. This should be compared with the sharper bound 
yWhex < 1.868832 obtained in [1] using the method of [2j with n = 45 and m = 17. 

It is clear that sharper upper bounds can be obtained for square and triangular 
lattices by Aim's method if carried out to higher values of n and m. Judging from 
the computational resources (928 CPU hours) required to obtain the bounds in |I. this 
should not be a very demanding calculation (compare this to the 25000 CPU hours used 
for the enumeration of the hexagonal lattice SAWs). 

4. Summary 

We have used Kesten's method of irreducible bridges to obtain improved lower bounds 
on the connective constant for self-avoiding walks on several planar lattices. The number 
of irreducible bridges is obtained by enumerating exactly the number of bridges using 
transfer-matrix techniques. In one approach we calculate the number of bridges of 
limited span but to great lengths while in a second approach we obtain an exact series 
expansion for the number of irreducible bridges. The first approach turns out to yield 
the sharpest lower bounds. The second approach allows us to study the convergence 
of the lower bounds fi(n) to the connective constant as a function of the truncation 
order n. We find that the limit is approached linearly in 1/n. In addition we use recent 
exact data for the number of SAWs c n to obtain some upper bounds on the connective 
constant. The upper bounds are generally much poorer than the lower bounds and also 
worse than those already obtained by other methods. We have summarised the results 
in Table H 
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